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Abstract 

We introduce a new family of orthogonal polynomials on the disk that has emerged in the 
context of wave propagation in layered media. Unlike known examples, the polynomials are 
orthogonal with respect to a measure all of whose even moments are infinite. 


1 Introduction 

For each a > — 1 there is a corresponding family of disk polynomials that are orthogonal with respect 
to the measure {1 — — y^)°‘dxdy on the unit disk D; these are sometimes referred to as generalized 

Zernike polynomials, named for the case a = 0 introduced in [8]. The well-established theory of 
disk polynomials is detailed in [SEIII]- The constraint a > — 1 stems from the requirement that 
the measure (1 — — y^)°‘dxdy have finite moments, which is necessary for meaningful evaluation 

of the corresponding scalar product 

{p,Q)a= / p{x,y)q{x,y){l-x^-y‘^Ydxdy (1.1) 

on arbitrary polynomials p and q. Recent work on the propagation of waves in layered media mm 
has brought to light a family of polynomials orthogonal with respect to (1 — — y'^)~^dxdy. Since 

f x^™'y^”^(l — x^ — y‘^)~^dxdy = oo (1.2) 

Jd 

for every pair of nonnegative integers m and n, the scalar product dni) is not defined for ar¬ 
bitrary polynomials in the case a = —1. Nevertheless, polynomials—which we term scattering 
polynomials—comprise an orthogonal basis for {p, dxdy/{1 — x^ — y^)). The purpose of the 
present paper is to present the details of this result. 


2 Definition and properties of scattering polynomials 


Referring to the notation z = x + iy for points in the unit disk D, one has the option of working with 
euclidean x, y-coordinates, or with complex coordinates z and 2 . As far as orthogonal polynomials 
are concerned these are essentially equivalent, as elaborated in [7]; the present paper uses whichever 
coordinates are most convenient for the task at hand. 

We define scattering polynomials by a Rodrigues type formula, as follows. For every (p, q) E 
with min{p, g} > 1, set 




(p,-?)(^) = 


_ ^ _(1 _ zz) (1 - zz)^+'?“b 

q{p + q — l)\ dzPdzi 


( 2 . 1 ) 
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The chosen normalization simplihes the formulation of the boundary Green’s function for scat¬ 
tering in layered media (see m) and so is physically natural, although not important for present 
considerations. 

Note that disk polynomials satisfy a Rodrigues formula similar to that of scattering polynomials, 
but there is a qualitative difference: it follows directly from (I2.1|) that = 0 for every 2 ; on 

the unit circle T, whereas all disk polynomials have constant non-zero modulus on T, cf. [3]. Our 
main result concerns completeness of scattering polynomials, as follows. 

Theorem 1 Scattering polynomials defined by 112.1\) . where {p,q) G and min{p, > 1, 

comprise an orthogonal basis for L‘^(D,dxdy/{l — — y^)). 

In ^2.11 and ^2.21 below we show that scattering polynomials are eigenfunctions of a second order 
differential operator and may be expressed in terms of Jacobi polynomials; these results contribute 
to a proof of Theorem [T] completed in ^2.,'ll 


2.1 Eigenfunctions of —(1 — — y^)A/4 

Let A denote the modified laplacian 


A = 


(1 - 2:z) 


dzdz 


1 — 


A, 


( 2 . 2 ) 


where A is the usual (euclidean) laplacian. Direct computation using (|2.ip shows that for all 
integers p,q > 1, 

-Aip^P’‘i')=pqip(P’‘i), (2.3) 

Letting uo : Z+ ^ Z+ denote the divisor function, there is thus a family of (To(A:) eigenfunctions 
of —A of the form corresponding to each positive integer eigenvalue k. We show in the next 

section that these eigenfunctions are linearly independent. 


2.2 Representation in terms of Jacobi polynomials 

Like disk polynomials, scattering polynomials have a representation in terms of Jacobi polynomials, 
but again, there is a qualitative difference. The disk polynomials corresponding to parameter 
a > — 1 can be expressed in terms of Jacobi polynomials for nonnegative integer values of 

fi. Since there is no Jacobi polynomial corresponding to a = —1, the same cannot be true for 
scattering polynomials. Indeed it turns out that scattering polynomials can be formulated in terms 
of Pn'^\ where /3 is a nonnegative integer, as follows. 

Expanding the binomial (I — zz)p^'^~^ in the formula (EU), and then applying the derivative 
dP+^ifdzPdzfi yields 






iz) = 




(1 - zz)z'^+’'-P+^z'^+’'-'i+^ 


j=0 


{j + n + m + I)! 
j\(j + m)\{u -j)\ 


(zzy, (2.4) 


where m = \p — q\ and n = min{p, g} — 1; the latter notation will be used in the remainder of this 
section. Switching to polar form 2 ; = re*®, it follows from (12.4h that 

ipM[re^0) =e*('?-P)®/(P’9)(r), (2.5) 
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where 


fM^r) = Yi-iy + ^ + + , 


. 2 i 


j=0 


( 2 . 6 ) 


The radial functions were hrst discovered in [2], as was the following connection to Jacobi 

polynomials, valid for p > 0 : 

fM(r) = + _ ^ 2 ymp{l,m)^ 2 r^ _ (2.7) 

Combined with (|2.5p this yields the representation 




Note that the angular part of (/?(P’'?)(re*®), namely 6 *^”^ is a pure frequency. Therefore if 
q — p y q' — p', then and are orthogonal in 


LCD 


dxdy 


1 — — y2 


= LG D 


rdrdO 
1 — 


In particular, if pq = p'q' and (p, g) / {p',q'), then and are orthogonal, so the set of 

scattering polynomials corresponding to any hxed eigenvalue of —A is linearly independent. 

2.3 Completeness in L^(d, Yr^ ) 

In general, given a measure // on a locally compact metric space X and a positive measnrable 
weight function tc : A —)• [R+, 

(2.9) 


I 


L^{X,wdp) = —=L‘^{X,dp), 


w 


and a sequence is an orthogonal basis for L?‘{X,w dp) if and only if the corresponding 

seqnence {^/wb^}'^Q is an orthogonal basis for Lp‘{X, dp). In particular, setting dp = rdrd6/{l—r‘^), 


L2(D, dp) = Vr^^L2(D, rdrdO). (2.10) 

Also, since for any nonnegative integer m, I Pi2’'^\u)\ is an orthogonal basis for 

I J u=0 

L‘^{[-l,l],{l-u)il + u)^du), 
it follows that the quasipolynomials 

T -u\^ (l + 


Q^Y\n) = 


comprise an orthogonal basis for L^([—1, l],(in); see [5]. 
In order to show that 


pa,m)( 


U 


( 2 . 11 ) 


13 = 


|(^(p>'?) I (p^gr) G & min{p, ( 7 } > l| 


( 2 . 12 ) 


December 24, 2014 














Peter C. Gibson 


Orthogonal polynomials on the disk 


4 


is an orthogonal basis of (D, rdrdO /(1—r^)), we first argue that the functions are orthogonal, 
and then that the span of B is dense. It was proven in ^2.21 that and p)^P'’^''> are orthogonal 

if pq = p'q' and {p,q) 7 ^ {p',q')- On the other hand, if pq 7 ^ p'q', then orthogonality of and 

i^ip tQ ) follows from the fact that they are eigenfunctions, corresponding to distinct eigenvalues, of 
the self-adjoint operator —A; self-adjointness of —A follows from that of —A by (I2.2p . 

It remains to show that span.^ is dense in L^(D, rdrd0/(l — r^)). Toward this end, suppose 
that h E L‘^{p,rdrd9/{I — r^)) is orthogonal to every member of B. By (|2.10p there exists g E 
L^(D,rdrd 0 ) such that 

h{r,6) = \/l — g{r,6). (2-13) 

Let ap,g = max.{p,q}/q denote the coefficient occurring on the right-hand side of 

(j2.8p . Then for each fixed n = q — p ^ Z, for every p = min{p, g} — 1 > 0, 

0 = [ Mr, 

= ap,, h{r, dO^ (1 - - 1 ) ( by m) 

= ap,q g{r, dO^ v^l - r 2 (2r^ _ 1 ) rdr (by (imiB 

= ^ iu = 2r^-l) 

= ^ ( as in dlTl])). 

Since the quasipolynomials are an orthogonal basis for L2([—1,1], drt), it follows that for 

each n G Z, 

[ = O (2.14) 

^ 0 

for every u E [—1,1] outside a set En of measure zero. Since is an orthogonal basis of 

L^([0, 27r], d0), it follows in turn that for u 0 ULln 


9 



= 0 , 


for almost every 9 E [0, 27r]. Thus g{r,9) = 0 for almost every {r,9) E D and g' = 0 as a func¬ 
tion in L2(D, rdrd^), whence h = 0 also. This proves that the orthogonal complement of B in 
L'^(J),rdrd9/{I — r^)) is empty, and hence that B is an orthogonal basis. 


3 Conclusions 

Since the vector space L‘^[D,rdrd9/{l — = y/l — r'^L"^ (D, rdrd9) is dense in 

L‘^[D,rdrd9) = (D, dxdy), 

and convergence in the former space implies convergence in the latter, scattering polynomials com¬ 
prise a (non-orthogonal) basis forL2(D, dxdy) consistent with Dirichlet boundary values. Prom the 
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perspective of analysis of functions on the disk, this provides an alternative to Zernike polynomials— 
and their generalizations the disk polynomials—which are non-zero on the boundary circle and so 
inconsistent with Dirichlet conditions. 

More generally, scattering polynomials illustrate that orthogonal polynomials can comprise an 
orthogonal basis for a function spaces in which not all polynomials are integrable. A 

natural question for further investigation is the existence and extent of other such examples. 
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